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O Abstract. This paper presents a parameter- free theory of shear-generated 

^^ turbulence at asymptotically high Reynolds numbers in incompressible fluids. 

^ It is based on a two-fluids concept. Both components are materially identical 

^ and inviscid. The first component is an ensemble of quasi-rigid dipole-vortex 

'^ tubes (vortex filaments, excitations) as quasiparticles in chaotic motion. The 

^ second is a superfiuid performing evasive motions between the tubes. The 
local dipole motions follow Helmholtz' law. The vortex radii scale with the 

' ^' energy-containing length scale. Collisions between quasiparticles lead either 

^ to annihilation (likewise rotation, turbulent dissipation) or to scattering 

Y (counterrotation, turbulent diffusion). There are analogies with birth and 

^ death processes of population dynamics and their master equations and 

^ with Landau's two-filuid theory of liequid Helium. For free homogeneous 

^ decay the theory predicts the TKE to follow t~^. With an adiabatic wall 

•'-H condition it predicts the logarithmic law with von Karman's constant as 

^ \I^2t{ = 0.399. Likewise rotating couples form dissipative patches almost 

i-^ at rest (— )■ intermittency) wherein under local quasi-steady conditions the 

i_i spectrum evolves into an "Apollonian gear" as discussed first by Herrmann 

^ [1990]. Dissipation happens exclusively at scale zero and at finite scales this 

j> system is frictionless and reminds of Prigogine's (1947) law of minimum (here: 

r^l zero) entropy production. The theory predicts further the prefactor of the 

Q 3D-wavenumber spectrum (a Kolmogorov constant) as |(47r)^/^ = 1.802, well 

lO within the scatter range of observational, experimental and DNS results. 

Q Keywords: Turbulence, vortex dipoles, vortex tubes, dipole chaos, two-fluid 

CN theory, quasiparticles, von Karman's constant, law of the wall, Kolmogorov 

T^ constant 



The diversity of problems in turbulence should not obscure 

the fact that the heart of the subject belongs to physics. 

- Falkowski and Sreenivasan (2006). 

1. Introduction the higher expansion elements subject to various clo- 
sure hypotheses [for a review see Wilcox, 2006]. This 

Many efforts to solve the turbulence problem rest on assumption is also reflected in one of the MiUenium- 

the idea that the Navier-Stokes equation (NSE) plays p^^e problems of the Clay Mathematical Institute an- 

the role of a God equation and the application of a nounced in 2000. However, until today these efforts 

certain number of mathematical operations onto NSE ^^^^d not answer most elementary questions about tur- 
could do it. In particular, the Fridman-KcUer [1924] 
series expansion of NSE played a prominent role, its 
zero*'' element being the Reynolds [1895] equation and 



bulence. Here we explain why^. 

A methodical alternative was chosen by Prandtl 
[1926] who discussed turbulence in terms of analo- 
gies with molecular diffusion, gas kinetics, and Brow- 
nian motion in the interpretation by Einstein [1905]. 
Prandtl related his mixing length (Mischungsweg) with 
the mean-free path of kinetic gas theory. This concept 
became popular but detailed questions could not be an- 
swered without use of measurements. Although Prandtl 
has been heavily criticized by Batchelor [1953], also 
other scientists derived free-hand analogies for turbu- 
lence, e.g. the early K-uj model by Kolmogorov [1942], 
corrected and improved by Saffman [1970], further im- 
proved by Wilcox [2006] . They are today part of a larger 
set of so-called two-equation turbulence models like K- 
e, K-KC, K-T etc.2. 

Three years before Prandtl, Debye and Hiickel [1923] 
based their theory of electrolytes on the assumption 
that each ion is surrounded by a spherical "cloud" (or 
screen) of ions of opposite charge. 

The ideas of Prandtl and Debye and Hiickel were 
early examples of constructive theories for large many- 
particle ensembles in very different branches of physics, 
where knowledge of parts and pieces does not suffice to 
describe the system's global behavior. 

A similar methodical challenge was met later in the 
fields of superconductivity and superfluidity by groups 
around Landau and Feynman [see Landau and Lifschitz, 
1983; Landau and Lifshitz, 1987; Feynman, 1955; An- 
derson, 2011]. They realized that their fundamental 
equations and principles^ hold perfectly for the single 
atom, but properties like electric conductivity (or the 
gloss of metallic surfaces) emerge only when a larger 
ensemble of atoms is put together so that new spe- 
cific features and interrelations come into play, without 
violating the fundamental equations and principles, of 
course. 

For superfluidity it was necessary to supplement (or 
even replace) fundamental equations by the two-fluid 
theory of Landau (1941) - a not too far relative of the 
ion-cloud concept of Debye and Hiickel [1923]. 

In many cases the formation of ensembles means to 
break symmetries. So eventually the concept of emer- 
gence and broken symmetries gave physics in the last 
years a new and much less rcductionistic face than it 



^Everywhere the pronomen "we" is used in this text, it means 
the two of us, the dear reader and the author. 

•^Comprehensive overviews on the history of turbulence theory 
can be found in Davidson [2004] and Davidson ct al. [2011] 

''Schroderinger equation, Pauli principle etc. 



had until 1923 [Glansdorff and Prigogine, 1971; Haken, 
1978; Laughlin, 2005; Anderson, 2011]. 

Encouraged by the victory of the a. m. two-fluid con- 
cept in superfluidity, Liepmann [1961; loc. cit. Spiegel, 
1972] and later on Spiegel [1972] made steps to test it 
for classical fluid turbulence. Along this road, which 
we follow here too, Spiegel [1972] was the very first 
who used terms like excitations, quasiparticles and two 
fluids in the context of classical fluid turbulence, fol- 
lowed later by Spalding [1985]^. These authors were 
at least fully aware of the principle, deep-rooted and 
almost philosophical problems with reductionistically 
NSE-bascd approaches. Ed Spiegel underlined his thor- 
ough conviction again just recently in the closing re- 
marks of the report on the final session of the Turbulence 
Colloquium Marseilles 2011^. Besides ancient precur- 
sors of our ideas like Rene Descartes (1596 - 1650), 
who spoke about "tourbillons" forming the universe, 
and Lord Kelvin, who coined 1867 the notion "vor- 
tex atoms" [loc. cit. Saffman, 1992] we have to men- 
tion IVIarmanis [1998] who possibly was the first to pro- 
pose vortex dipoles as the fundamental quasiparticles of 
turbulence. Finally, the numerical JVIonte-Carlo eddy- 
collision methods with various vortex-filament primi- 
tives indicate that ideas developed by practitioners are 
not too far away from our theoretical views [e.g. An- 
deme, 2008]. 

Below we elaborate the two-fluid concept in greater 
detail; not exhaustive because the number of potential 
applications and side-problems is huge. But we show 
that, in an idealised sense, turbulence in an incompress- 
ible fluid at i?e — >■ oo can be understood as a statistical 
many-body ensemble - a tangle of vortex-dipole tubes 
(or filaments) taken as interacting discrete particles. We 
will answer a number of open questions of turbulence 
without use of empirical parameters. 

2. Broken symmetry and irreversibility 

In his anti-reductionistic article summarizing experi- 
ences of condensed-matter physics, Anderson [1972] 
states that "more is different". The turbulence theory 
presented below is an extreme example. While the case 
of one quasiparticle (an almost frictionless vortex dipole 
of zero net circulation) in a volume is non-dissipative 
and symmetric with respect to time and circulation, al- 
ready the presence of two quasiparticles in a common 
volume has the potential to break the symmetry and 



^see www.archive.org/details/nasa/_techdoc/_19860008217 
^See http : //turbulence . ens . f r/ 



to give rise to the emergence of turbulent dissipation: 
according to the laws of Helmholtz, dipoles are always 
in motion, may thus collide and may - depending on 
the occasional collision angle - form likewise-rotating 
and thus unstable couples. Their centers of mass stay 
almost at rest and evolve into a spectrum of smaller and 
smaller vortices until their kinetic energy is converted 
into heat at scales of size zero. In a sense we have here 
a most simple many-body problem because already the 
transition from M = \ to N = 2 suffices to break the 
symmetry and allow the emergence of dissipation. 

Setting Re — > oo in NSE means vanishing viscosity 
such that only the Euler equation remains. However, 
the latter has non-unique solutions, visible already in 
the trivial case of the inviscid Burgers equation. Ad- 
ditional information is needed to achieve uniqueness. 
I.e. with Re — > oo we first delete physics (viscosity) in 
NSE to be forced then to add (from outside) reasonable 
physics to finally reinstate uniqueness. But this is not 
all to be overcome. In addition, we have the problem 
of localization and integrability of a solution. Classi- 
cal weak vortex solutions of the Euler equation extend 
into the full volume and are thus not integrable. But 
real-world vortex ensembles exhibit individually local- 
ized vortices and finite scales. 

Those contradictions could only be resolved intro- 
ducing quasiparticles embedded in the two-fluid concept 
initiated earlier by Landau, Tisza, Feynman, Liepmann, 
Spiegel and Spalding, resting all on Debye and Hiickel 
[1923]. 

3. Vortex tubes and dipoles: Batchelor 
couples and von-Karman couples 

Potential vortex. The classical potential vortex 
around a closed vortex line (e.g. the centerline of a 
smoke ring) represents an exact weak solution of the 
Euler equation. The vortex line has infinitely thin di- 
ameter, infinitely high angular velocity, but finite cir- 
culation, c = TT p^ oj < oo, where p and cj are radial 
coordinate and vorticity. The fluid outside the vortex 
line is inviscid and irrotational and the radial velocity 
V — LJ p ^ 1/ p decreases with distance p from the cen- 
terline and extends into space. It is not applicable here, 
as well as his relative, the Rankine vortex, because they 
both are filling the space. 

Vortex tube. Instead of the above we use the 
vortex-tube concept, which is a Rankine vortex without 
sticking condition between its forced and free parts. I.e. 
its radial velocity increases (like in a rigid vortex) lin- 
early from center to tube radius r, and for p > r there 



is an inviscid and independent potential fiow governed 
by volume conservation. Thus vorticity is confined to 
the interior of a spaghetti-like tube as described e. g. by 
Lugt [1979] and in greater detail by PuUin and Saffman 
[1998] who quote papers by Kuo & Corrsin, 1972, and 
Brown & Roshko, 1974, showing tubes as dominat- 
ing characteristic structures. A more recent simulation 
study has been presented by Wilczek [2011] in form of 
a number of instructive movies on vortex ensembles in 
motion on the internet®. 

The centerlines of our vortex tubes form either closed 
loops or they are attached to boundaries. The problem 
of stability of the tubes is discussed further below. 

Batchelor couple. This is a vortex dipole made up 
of two anti-parallel vortices [Lesieur, 1997] and symbol- 
ized by (+1)" — ) or (— f|- +) where plus and minus mean 
the signs of vorticity within the vortices and the arrow 
the direction of motion of the couple or dipole. In clas- 
sical interpretation the flow field of one vortex moves 
the other vortex and vice versa. The total circulation 
of a dipole is zero because the vorticities carry opposite 
signs. In practice such couples are stable over mod- 
erate propagation times. In our idealized image they 
are made up of vortex tubes, move frictionless with lo- 
cal center-of-gravity velocity u — Lur and conserve all 
their properties excepting their position in space be- 
cause they propagate with their local u into the same 
direction as the fluid between the two tubes. In a dense 
ensemble of Batchelor couples their trajectories are no 
longer straight lines due to mutual interactions. A cou- 
ple's kinetic energy density is m^/2 = r^u;^/2. 

Von-Karman couple. It is a counterpart of the 
Batchelor couple, made up of likewise rotating vortex 
tubes and symbolized by (+ || +) or (— || — ). Its to- 
tal circulation does not vanish. Such a couple is known 
since long to be fundamentally unstable. Its kinetic en- 
ergy is eventually dissipated into heat [e.g. Lamb, 1932]. 
Further below we discuss details of this process. 

4. Dipole chemistry in reaction-diffusion 
approximation 

For 2D trajectories it has been found by Aref [1983] 
and Eckhardt and Aref [1988] that the trajectories are 
chaotic so that for very high Re and 3D motions of 
tube-like vortex dipoles in form of a dense 3D tangle 
we may assume also chaotic motions where collisions 
cannot be excluded. Let us consider an asymptotically 

^^http: //pauli .uni-muenster .de/tp/menu/f orschen/ 
ag-f riedrich/mitarbeiter/wilczek-michael . html 



high tube density so that the chaotic trajectories be- 
tween colhsions are short and locally homogeneous (in 
a statistical sense), like e. g. in the theory of Brownian 
motion. Then the total dynamic process of the tan- 
gle may be described as [Frisch, 1995] ". . . conservative 
dynamics punctuated by dissipative events. . . " : some 
"elastic" collisions lead to energy-conserving reorgani- 
zations of Batchelor couples, resembling turbulent diffu- 
sion, whereas other collisions are dissipative or "inelas- 
tic" and lead to the formation of fundamentally unsta- 
ble von-Karman couples whose energy moves to smaller 
and smaller scales where it decays, resembling turbulent 
dissipation or dipole annihilation. 

free di poles with 
zero total circulation: 
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recombirred stable dipoles: 
turbulent diffusion 



instable configurations: 
turbulent dissipation 



Figure 1. Collisions of two dipoles: the left path- 
way is "diffusive"; it is a recombination of dipole el- 
ements and chaotically scatters the trajectories (turbu- 
lent diffusion). The right pathway is "dissipative"; it 
evolves into an unstable vortex configuration which de- 
cays "somehow" into heat. 

Figure 1 shows the two potential results of a dipole- 
dipole collision. For symmetry reasons the two path- 
ways have identical probabilities of 72- 

The simplest mathematical structure describing the 
statistics of the above two irreversible processes of (tur- 
bulent) diffusion and (turbulent) dissipation of dipoles 
is given by the following reaction-diffusion equation, 
which may be understood also as a special case of the 
Oregonator [see e.g. Ch. 9 in Haken, 1978]: 



— — fuAf - — 
dt dx \ dx 



= F-PN^ 



(1) 



Here N is the volume density of dipole tubes (or fila- 
ments), V is turbulent diffusivity, /3 is a constant, F is 
a source term describing the generation of dipole tubes 
and set equal to zero for the moment, v and /3 are un- 
known so far. We only know for sure that the exponent 
of M in the last term of (1) is really 2 because it needs 
two colliding dipoles to generate (with probability Y2) 
unstable configurations which dissipate energy. But if 
energy is gone, a whole dipole is gone as our dipole tubes 
as quasiparticles differ from their inviscid potential-flow 
environment only by their kinetic energy. 

The presence of a mean flow and the corresponding 
advection of dipoles with the flow is sketched for reasons 
of completeness by the term with mean-flow vector U . 
For simplicity we set further always [/ = 0. 

One may look at (1) only from the viewpoint of pure 
analogy with chemical reaction-diffusion processes. But 
one may also derive this equation over many pages from 
scratch, beginning with a master equation for the prob- 
abilities as sketched earlier [Baumert, 2005b, Ch. 5.6]. 
Fortunately, this has already been done by other au- 
thors many years ago for whole classes of such pro- 
cesses, and went into the many textbooks on stochastic- 
dynamic systems, physical kinetics and other fields 
[e.g. Kraichnan, 1968; Haken, 1978, 1983; Stratonovich, 
1992, 1994]. Therefore we will not repeat here these 
formal steps and refer the more technically interested 
reader to the literature mentioned. 

5. Two fluids, dressed and naked tubes 

Two fluids. We consider a circulation- free volume 
filled with two different forms of an materially identical 
incompressible fluid: 

(a) inviscid, quasi-rigid and incompressible but de- 
formable vortex tubes as quasiparticles with finite 
radius r and vorticity module w, vorticity bundled 
within the tube; in our analogy with the kinetic 
theory of gases a dipole tube is an analogue of a 
particle, 

(b) the inviscid fluid between the dipole tubes; it is 
and in the analogous kinetic theory of gases an 
analog of the vacuum. 

Fluid (b) behaves like a super-fluid and receives no force 
from the moving quasi- rigid dipole tubes (d'Alambert's 
paradox). The fluids differ only in their state of mo- 
tion. While the tubes rotate around their (in general 
curvilinear) axes and move locally relative to the vol- 
ume according to Helmholtz' laws, the fluid between the 



tubes performs corresponding evasive motions accord- 
ing to the principle of volume conservation. 

Naked tubes. Above we have used the concept of 
vortex dipoles made up of ideal vortex tubes immersed 
in an inviscid fluid and exchanging no energy with it. 
However, this concept is only applicable if the vortex is a 
rigid body. If it is a fluid, the problem of stability arises 
because the quasi-rigid vortex tube as an exact solution 
of the Euler equation is accompanied by the following 
pressure head as a consequence of incrtial (centrifugal) 
forces: 

p = Po + ^xw^r\ (2) 



Here Pq is the background pressure of a laminar refer- 
ence flow, e.g. in the ocean the depth-depending hydro- 
static pressure. If the pressure outside the vortex would 
simply be po then, due to the action of the outwards- 
directed pressure force given by the second term in (2), 
the vortex would loose stability against small pertur- 
bations. However, tubes of finite radius are observed 
to be during their finite life time relatively stable in 
real- world turbulence. 

Dressed tubes. This contradiction can be ex- 
plained by the consideration of ensemble effects. It has 
been noted in the Introduction of this paper that in 
many-body problems like turbulence the phenomenon 
of emergence deserves special attention. This means 
that in a (local) volume element with a larger number 
{N ^ 1) of similar vortex tubes (a local "cloud") the 
tube ensemble itself generates the background pressure 
(2) which keeps eventually all the individual tubes - 
at least in the center of the cloud - in a "sufficiently 
stable" state. 

In thermodynamically open systems like turbulence 
such a cloud represents a quasi-steady state [FUefigle- 
ichgewicht in the sense of Bertalanffy, 1953; Glansdorff 
and Prigogine, 1971; Haken, 1978, 1983]. I.e. the pro- 
cesses of dipole generation (see below) and their annihi- 
lation by collisions almost compensate each other. Any 
quasiparticle (dipole) has thus only a limited statisti- 
cal life time. Therefore "sufficient stabilization" of a 
dressed dipole tube by a cloud means to guaranty sta- 
bility only in a statistical sense during its (statistical) 
life time. 

6. Particle number, TKE, 

and r.m.s. vorticity frequency 

Consider a small volume element 6V populated by an 
ensemble of j = 1 ... A/" dipoles with individual effective 
vortex radii, rj, and r.m.s. vorticity moduli, ujj. The 



latter can be interpreted as r.m.s. values of individual 
dipoles j as follows. 
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arc the individual vorticitics of the 



where +ujj and —ljj 

two vortex tubes forming the dipole. 

The dipole energy is conserved as long as dissipative 
events do not take place. The volume density of dipoles 
is Af/SV. The total TKE within 5V is the sum of the 
kinetic energies of the individual dipoles: 
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Multiplication of (4, 5) with 5V gives 
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K = SV-ICsv 
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Ksv, ^sv are local volume densities of TKE and r.m.s. 
vorticity magnitude, respectively. They and AT, w are 
extensive variables by definition, i.e. they scale with the 
dipole number inSV. k and w are ensemble averages and 
as such intensive variables which do not change when 
new particles with average properties are added to SV: 



k = 



-y - 



(jj 



Af ^ 2"^"^ 
jesv 

jesv 



(8) 
(9) 



We turn now to equation (1) where, according to (6, 7) 
we replace M with Kjk and w/w, respectively, to get 
eventually balance equations for the extensive variables 
K and fJ = a;/27r, provided that the intensive variables 
k and uj vary sufficiently weakly in time and space com- 
pared with the dipole number M: 



— —ivK-v — 
dt dx \ dx 

dt dx \ dx 






Here fl = l/T is the ordinary^ vorticity frequency and 



^While uj = 2tt/T is an angular frequency, f! 
ordinary frequency. 



l/T is an 



related with w by the often used constant k which fur- 
ther below appears to be von-Karman's constant: 



w = 2nn = n/n^, 


(12) 


K = (2 7r)-i/2« 0.399 . 


(13) 


see that 




f^K = P/'k, Pn = 27r/3/w, 


(14) 


Tk = kT, Tn = wJ"/27r. 


(15) 



7. Mixing length and eddy viscosity 

Prandtl [1925] hypothesized that his mixing length has 
a two-fold meaning. It should be "considered as the 
diameter of the masses of fluid moving as a whole" (in 
our picture: as a characteristic vortex radius) or "as 
the distance traversed by a mass of this type before it 
becomes blended in with neighboring masses. . . " (in 
our picture: as a mean free path), [loc. cit. Bradshaw, 
1974]. I.e. he stated a tight relation between the radius 
and the free path of a dipole in a dense ensemble. In- 
deed, for a dipole tangle this can be shown. We first 
define the ensemble average, f, of the vortex radii as an 
average over a volume element weighted by vorticity as 
follows, 
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(16) 



We now determine the eddy viscosity, v^ in analogy to 
the theory of Brownian motion [Einstein, 1905] in terms 
of a mean free path. A, and a mean free flight time, 
T = \ju: 



V = AV2t. 
Here u is the mean velocity of a dipole. 



(17) 



= 2/f = Vti? = Vr?c^?. (18) 
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For a turbulent tangle of dipole tubes at Re = cxd 
in a quasi-steady state far from solid boundaries the 
mean free path will clearly be short but cannot vanish 
because there is an equilibrium of dipole generation, 
annihilation, and motions that lead to the latter. This 
implies that [Baumert, 2005b] 

A = If. (19) 

Algebra finally gives the following: 

V = K/nVt. (20) 

This relation is called the Prandtl-Kolmogorov relation. 
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Figure 2. Local cross section through a dipole tan- 
gle far from boundaries at maximum dipole density 
{Re — >■ oo). The dark grey dipole in the left half square 
labels the initial situation, the right half square a situ- 
ation after a translation motion from left to right into 
its collision position, where it is either annihilated or 
scattered. A free dipole obviously occupies statistically 
a cross sectional area of (4f)^. 

8. Parameters /3x, /Jn, and C, 

We consider an initial-value problem {Fk — 0, Fn = 0) 
in a spatially homogeneous volume where dK/dx — 
and dfl/dx = such that equations (10, 11) reduce to 



dK 

~dt 

dn 

~dt 






0, 
0. 



For large t follows that 

Kit) = {PKt)-\ 
^(t) = {Pnt)-\ 

and, according to (20), for eddy viscosity holds 

K{t) Pn 



^{t) 



TTn{t) 



ttPk 



const. 



(21) 
(22) 



(23) 
(24) 



(25) 



Equation (23) coincides with the results of a fairly gen- 
eral similarity analyses of the Navier-Stokes equations 
by Oberlack [2002] and with experimental results by 
Dickey and Mellor [1980]. 

We now use the definition of the dissipation rate, e, 
which is the last term on the right-hand side of (10) 
and carries the units of of TKE (density) per time, i.e. 
e ^ K/T ^ Kil, [m^ s~'^]. For reasons of convenience 
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we write this variable as follows and introduce a still 
unknown auxiliary variable, C: 

t2 



PkK' 



(■ 



Kn 



C- 



(26) 



where we made use of (20). This implies 



(27) 



When the dipoles behave in general like in free decay 
then we may use (25) which gives then the expressions 

Pn = C/tt, (28) 

(3k = C/i^t:^. (29) 

With these results we may rewrite (10, 11) as follows: 
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dx J 
dfl 
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(30) 



Co2 



= Fn- -n\ (31) 



dx J IT 

with V given by (20) and (^ = 1 as shown further below. 

9. TKE and vorticity generation 
by mean-flow shear 

Generation of TKE. To further complete equations 
(30, 31), we have, besides ^, to specify the source terms 
^K , Fn ■ The specification of Fk for shear flows is triv- 
ial because it is given by the classical losses of the mean- 
flow energy balance and can therefore be copied from 
textbooks like Wilcox [2006] or Schlichting and Gerstcn 
[2000]: 



J"k = - ^ (u ■ Mj- 
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•'' dxi 



- -S'-'-Kj26. 



i,j = i 



dxi 



(32) 



where — (u^u') is the Reynolds stress tensor defined as 
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and 5*,, is the rate of strain tensor defined as 



bj. 



dUi dU, 
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dxi dxj 



S^ is the total instantaneous shear squared, 
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i,i = l 



dU, 
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(33) 



(34) 



(35) 



Ui is the z*'* component of the mean flow velocity vec- 
tor U = {Ui, U2, C/3)"^ and 6ij is the Kronecker symbol 
which is zero for i ^ j and unity for i = j. In the simple 
case of a plane horizontal flow with vertical shear like 
wind over plane terrain or flow in a plane wide channel 
we have i = 1 and j = 3 such that (32) reads as follows: 



Fk = i^S^ 



(36) 



Generation of vorticity. We use a fundamen- 
tal macroscopic argument by Tennekes [1989]. It has 
flrst been cast into mathematical form by Baumert and 
Peters [2000, 2004] and carefully discussed by Kantha 
[2004, 2005], by Kantha et al. [2005] and by Kantha 
and Clayson [2007]. Tennekes hypothesized that, in a 
neutrally stratified homogeneous shear flow, an energy- 
containing turbulent length scale, C, cannot, on dimen- 
sional grounds, depend on the ambient shear. 

Further, also on dimensional grounds we have L ~ 
K^/^n~^ and thus L^ ~ K^'"^ such that the time 
evolution of the length scale follows 



1 dL^ 



1 dK 
K lU 



1 dn 



(37) 



We replace dK/dt with Fk according to (32) and flnd 

(38) 



1 dL^ 



S^ 



1 dn 

2 

n dt 



Tennekes' hypothesis means that the evolution of i can- 
not be controlled by S which means dL^/dt — and 
implies that 

^ = ^. (39) 

where we replaced dVt/dt with Fu. Algebra gives 

572 TT (40) 



Fo. - 



and for the simple case of a plane horizontal flow with 
vertical shear we may complete with some algebra equa- 
tions (30, 31) as follows: 



= i^{S^- C02) , (41) 
1 /52 
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— + — un- 

dt dx V 
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- C^' (42) 



with 1/ again given by (20). 

10. Turbulent boundary layer 

Boundary conditions. Consider a horizontally homo- 
geneous flow and its stationary boundary layer close to 



a plane solid wall at X3 = z = where for convenience 
z is introduced here as the only relevant coordinate, 
pointing orthogonal from the wall into the fluid. Thus 
the shear is 

S = \dU/dz\ . (43) 

This special case means that in (73, 74) the horizontal 
and time derivatives vanish and it remains 
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dz 

d 
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dK 
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i^{s^-(:n^), (44) 

C^A (45) 
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The diffusive TKE flux into the viscous sublayer a.t z — 
Zq has to vanish in the sense of an adiabatic boundary 
condition, 

{lydK/dz)^^^^ = 0, (46) 

such that also the flux divergence on the left-hand side 
of (44) vanishes, giving K ^ K^ and v {S"^ - C, il^) = 
or 



n = s/^(. 



(47) 



Logarithmic law of the wall. We insert (47) into 
(45). The unknown C, cancels out^ and we have to solve 
the following equation for S = S{z), 



2K,l(l-.'^\ ^ S' 



dz \ S dz 



which gives 



S{z) 



dU V2 Ka 



(48) 



(49) 



dz z 

Integration of (49) gives the logarithmic law of the wall. 
In boundary layer theory the bottom shear stress is 

dU Ko 



defined in terms of the squared friction velocity, ui , 



•* dz TT O ' 



and with (47) it follows that 

This allows to rewrite (49) as follows, 

dU Uf 

dz kz' 



(50) 



(51) 



(52) 



*The respective text in http : //arxiv . org/pdf /0907 . 0223 . pdf 
contains an algebraic error but without later consequences, for- 
tunately. 



with K as a modified von-Karman constant defined with 
respect to (15) through 



Integration of (52) provides us with 



U(z) 



Uf 



In 



(53) 



(54) 



Mixing length L. Consider the definition of the 
effective (statistically averaged) dipole radius of an en- 
semble through (16). We solve this equation for K and 
express the TKE in terms of f and fl as follows: 



K 



2tt 



2 ^2 



n' 



(55) 



Following now Hinze [1959, p. 279, eq. 5-2], in present 
notation Prandtls mixing length L, which is termed also 
the "energy-containing length scale" in the literature, is 
defined in terms of eddy viscosity and shear as follows: 



L' 



dU 
dz 



L'S. 



(56) 



Due to the eddy-viscosity formula (20) relation (56) 
gives 

so that in the neighborhood of a solid wall we get with 
(47) the following result: 



Comparing (55) with (58) gives 
L = f/k. 



(58) 



(59) 



The physical meaning of L can be understood as follows 
(see Fig. 3). If we may set C = 1 then (59) gives together 
with (15) and (53) the following relation: 



L^ = 2x (7rf2). 



(60) 



It means that L is the length of a square with an area 
equal to the cross sectional area of a dipole (in a statisti- 
cally averaged sense) because w f ^ is the cross-sectional 
area of one vortex tube. In an asymptotic sense this case 
corresponds to the maximum deformation of a dipole 
and justifies to set C = 1- 

While Prandtl's mixing-length concept was applica- 
ble only in the vicinity of solid boundaries so that it at- 
tracted respectful criticism Wilcox [2006] our concept is 
a generalization of Prandtl's concept and works also far 



from boundaries, even in the free stream of stratified flu- 
ids where L may approach the Thorpe scale and/or the 
Ozmidov scale, depending on the conditions (Baumert 
and Peters 2004). 

We summarize this section as follows: 



L 

z 



UfL, 

K Z , 

L/k = f/n^ . 



(61) 
(62) 
(63) 




///\\< 



VTT^S 



Figure 3. Cross section through a dipole sheet at a 
solid boundary. In this Figure we have r = f, i.e. the 
sketched quasiparticle is to be understood as an ensem- 
ble average. 



11. Spectra, dissipative patches, 
and spectral constants 

Spectra. Up to now our discussion was concentrated 
on only a few scales like f and L interrelated by n, and 
on T interrelated with Q and uj etc. However, the reality 
of turbulence exhibits a whole range of scales at which 
fluctuations occur. They have relatively stable spectral 
properties. This problem has attracted early attention 
by Richardson [1922] and Taylor [1937]. Kolmogorov 
[1941] found on dimensional grounds that the kinetic 
energy spectrum as function of wavenumber, wherein 
energy flows steadily from the large (energy-containing) 
scales to the much smaller dissipative scales, may be 
written as follows. 



Here k = 27r/A is the wave number and A the wave- 
length, e is the dissipation rate of TKE. Dimensional 
arguments force that that a2 = 2/3 and as = 5/3, 
in agreement with the oceanographic observations by 
Grant et al. [1959] in a tidal inlet with Re « 10* and a 
depth of about 100 m. A value for ai is derived below. 
Devil's gear. Our view of the above Kolmogorov- 
Richardson cascade has been filled with life through a 
numerical simulation study by Herrmann [1990] who 
demonstrated that Kolmogorov's value for a^ corre- 
sponds numerically to the data of a space-filling bearing 
[see also Herrmann et al., 1990]. The latter is the dens- 
est non-overlapping (Apollonian) circle packing in the 
plane, with side condition that the circles are pointwise 
in contact but able to rotate freely, without friction or 
slipping [a "devil's gear" sensu Poppe, 2004]. The con- 
tact condition for two different "wheels" with indizes 1 
and 2 of the gear reads 



U = UJiri = UJ2 ^2 



(65) 



where u is necessarily constant throughout the gear and 
governed by the energy of the decaying (initially energy- 
containing) vortex pair as u = \/lK . It follows that 



W2 = Wi 



^2 



(66) 



dK = ai e"^ /c""^ dk . 



(64) 



and for very small r2 the frequency a;2 may become 
high, even acoustically relevant. 

The dissipative patch. If the above gear is fric- 
tionless then the question arises where energy can be 
dissipated. In a real fluid with non-zero viscosity, dissi- 
pation happens at all scales, mainly but not exclusively 
where the velocity gradient is highest, here: at a scale 
vanishing with Re — >■ cx) to the size zero. Our dissipa- 
tive patch (Fig. 4 shows the first stage of its formation) 
is thus "almost frictionless" and a Hamiltonian clock- 
work, excepting scales of size zero. 

The formation of a fully developed spectrum of 
"wheels" from Fig. 4 deserves certain perturbations 
"from the sides" , a condition which is guaranteed by the 
random reconnection/recombination and scatter pro- 
cesses sketched in the left half of Fig. 1 and also by the 
incomplete mutual pressure compensation of the vor- 
tices in our vortex ensemble. In a quasi-steady state 
these perturbations initiate roll-up instabilities at the 
boundaries of the respective larger vortices so that even- 
tually and in a statistical sense a patch like in Fig. 4 is 
formed and evolves steadily into a fully developed gear. 

Our Fig. 1 illustrates the possible results of a dipole- 
dipole collision. While the left pathway shows the re- 
combination of coMTiter-rotating vortices from counter- 
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rotating vortices, the right shows the formation of a cou- 
ple of likewise rotating vortices from counter-rotating 
vortices. The latter then revolve around a common 
center of mass which remains nearly at rest (Fig. 4). 
Such a couple is fundamentally unstable [Sommerfeld, 
1948], a quasi-steady dissipative patch evolving into a 
full gear in the sense of the mechanisms discussed in the 
last Chapter. This picture lets us expect that dissipa- 
tion should exhibit a spatially patchy behaviour which 
we may also call intermittency. This problem has been 
studied extensively by various authors from other points 
of view [see e.g. Frisch, 1995] and cannot yet be dis- 
cussed here from our viewpoint in greater detail. 

Kolmogorov's constant cti. This constant be- 
longs to the wavenumber spectrum and deserves an idea 
about the outer limits and inner structure of an unsta- 
ble, dissipative patch as sketched in Fig. 5 for the begin 
of the cascade process evolving into a strucure like the 
one given in Fig. 4. 

The most important message of Fig. 5 is that the 
longest or energy-containing wavelength of the dissi- 
pative patch equals Ag = 2f. The wavelength in a 
dipole is 4f. The dipole performs chaotic trajectories 
in a white-noise sense and forms no patch or spectrum. 
This difference between the two configurations is essen- 
tial. We use our Aq as a lower integration limit for 
the spectral energy distribution. It is important to un- 
derline that Aq labels the upper wavelength limit (the 
longest wavelike motions) in a dissipative patch. This 
limit is actually not influenced by the formation details 
of the spectrum. 

We integrate (64) over the dissipatve patch in the 
sense sketched in Figures 4 and 5 and get 



K = aie 



2/3 



i>-'"*-'lfe)"' 



(67) 



where ko — 2 tt/Ao characterizes the lower end of the 
turbulence spectrum in the wavenumber space. We 
loosely assign the wavenumber range k = . . . ko to the 
mean flow which may basically be resolved in numerical 
models. The dissipation rate e in (67) can be expressed 
as follows, 

£ = K/t , (68) 

with r being the lifetime of a dissipative patch. Insert- 
ing (68) in (67) and rearranging gives the following: 



.,.|,2,)-A->«(^) 



2/3 



(69) 




Figure 4. Cross section through the first developmen- 
tal stage of a dissipative patch, i.e. of an unstable pair of 
likewise rotating vortices resulting from a dipole-dipole 
collision (right pathway in Fig. 1). The green circles 
represent primary energy-containing vortices. They do 
not touch each other due to spontaneously formed sec- 
ondary vortices which initiatiate a whole vortex cas- 
cade. The broken blue line and the arrows symbolize 
the slow rotation of the patch around its center of mass. 



K. 


TY' /I 
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Figure 5. Outer limits of a dissipative patch (c./. Fig. 
4). The maximum wavelength is obviously equal to 
Ao = 2f. 



places which are simultaneously also "emptied" through 
dipole annihilation or dissipative patches by decay. This 
means that the life time of a dissipative patch, t — K/e, 
should equal the time of "free flight" of a dipole over a 
distance 2 f: 

T = K/e = 2 f/u . (70) 



In a local quasi-equilibrium sense for a dense vortex 
ensemble the marching dipoles can occupy only those 
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Here we used the scalar dipole velocity u, 
u = u! r = V 2 K. 



(71) 



After some algebra we get the pre-factor of the three- 
dimensional wavenumber spectrum as follows: 



ai 



-(47r)2/3 = 1.802. 



(72) 



The corresponding value of an ideal one-dimensional 
spectrum is one third of the above, i.e. 0.60. 

12. Discussion 

Equations of turbulent motion. The results of our 
considerations can be summarized as follows: 



dK 

1h 



d_ 
dx 
d_ 
dx 
K 



U K -V 



un- 



dK 

dx 

on 

dx 



= iy{s^- n^ 



1 fS' 



-n^ 



(73) 

,(74) 
(75) 



These equations are structurally identical with the k-uj 
closure model discussed by Wilcox [2006]. There are 
only slight differences in the pre-factors of the terms. 

This theory applies exclusively to locally homo- 
geneous, isotropic and moderately unsteady high- 
Reynolds number flows. Extreme non-stationarities 
and/or sharp spatial gradients like in Shockwaves are 
possibly not covered. As a rule, temporal changes of 
the mean flow should happen on time scales sufficiently 
long compared with T — 1/17 because otherwise spec- 
tral universality (64) has possibly not enough time to 
become well enough established. 

Our equations reproduce the logarithmic law of the 
wall^ and predict the universal von-Karman's constant 
as K = l/-\/27r = 0.399 « 0.4 where the latter value 
counts as internationally accepted standard [see also 
Hogstrom, 1985]. The value corresponds nicely to mea- 
surements under favorable pressure gradients [Chauhan 
et al., 2005]. Similar support comes from Jimenez and 
Moser [2007] on the basis of an extensive review. They 
state: The Kdrmdn constant n Ri 0.4 is approximately 
universal. 

However, Landau and Lifshitz [1987] wrote on their 
p. 173 that . . . K (is) a numerical constant, the von- 
Kdrmdn constant, whose value cannot he calculated the- 
oretically and must he determined experimentally. It is 
found to be K ~ 0.4. 



^Thc controversy of the logarithmic versus a power-law 
boundary layer [Barenblatt, 1997] shall not be discussed here. 



With respect to the measurability of k even the op- 
posite might be true: if "physics disappears" when 
Re — > oo and only sort of "inert geometry" (the Eu- 
ler equation) remains, k can possibly no longer be seen 
as a physical quantity as it represents pure geometry. 

But these interesting speculations shall not rise any 
doubt about the necessity of superpipe facilities for very 
high-i?e number experiments like in Princeton [Zagarola 
and Smits, 1998], Oregon [Donelly, 1998] or in the Eu- 
ropean CICLoPE [Talamelh et al., 2009]. 

Our equations describe the free decay of turbulence 
following K ^ i^™ with to = 1, in agreement with 
Dickey and Mellor's [1980] high-i?e laboratory experi- 
ments and with Oberlack's [2002] theoretical result for 
the Navier-Stokes equation. Today it is still not clear 
why some decay experiments lead to to > 1. Possi- 
bly it is a matter of initial conditions [see Hurst and 
Vassilicos, 2007]: at high Re viscosity is comparatively 
small so that its regularizing effect towards a finally 
more self-similar decay spectrum will take more time 
than at lower Re. In some cases this time may exceed 
the lifetime of turbulence. This underlines the necessity 
of deeper experimental work. 

Kolmogorov's constants. The rounded numerical 
values ai = 1.8 or ai/S = 0.6 predicted by our the- 
ory for Kolmogorov's universal constant are situated 
well within the error bars of many high-i?e observa- 
tions, NSE and RG based analytical approximations, 
laboratory and DNS experiments. 

Based on observations, Tennekes and Lumley [1972] 
gave the value ai — 1.62, but still with some uncer- 
tainty. The study by Sreenivasan [1995] (see Fig. 6) is 
possibly the most comprehensive review of experimen- 
tal and observational values for the number ai/3 until 
now. Later Yeung and Zhou [1997] reported a value of 
ai = 1.62 based on high-resolution DNS studies with 
up to 512'^ grid points. A most recent study by Donzis 
and Sreenivasan [2010] on a DNS grid of 4096^ gave 
ai « 1.58. 

Much higher Reynolds numbers than in DNS could be 
achieved in oceanic measurements of Lagrangian fre- 
quency spectra by Lien and D'Asaro [2002]. These au- 
thors stated for the prefactor /3i in the Lagrangian fre- 
quency spectrum that . . . since the present uncertainty 
is comparable to that between high quality estimates of 
the Eulerian one- dimensional longitudinal Kolmogorov 
constant measured by many dozen investigators over the 
last 50 years, large improvements in the accuracy of the 
estimate of /?i seem unlikely. For completeness reasons 
we mention other theoretical efforts to calculate the uni- 
versal constants. They are technically extremely com- 
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Figure 6. Experimental and observational results 
for ai/3 measured, collected from the literature, and 
analyzed by Sreenivasan [1995]. The solid line fol- 
lows our somewhat arbitrary approximation 0.6 x 
y/Re\/[y/Re^ + y/Re\) wherein ai/3 = 0.6 is our the- 
oretically derived asymptotic value. Here we took 
i?e* = 10. 



plex and neither unique nor part of an integrated de- 
scriptive concept for turbulence so that they all carry 
more singular characters. Beginning with an initiating 
work by Forster et al. [1977], systematic analytical ap- 
proximations using RG methods and related techniques 
for NSE gave rise to some estimates. E.g. Yakhot and 
Orszag [1986b, a] found ai « 1.62 whereas McComb 
and Watt [1992] derived ai = 1.60 ±0.01 and Park and 
Deem [2003] obtained ai = 1.68. 

Coda. Saffman [1977] [loc. cit. Davidson, 2004, p. 
107] feared that . . . in searching for a theory of turbu- 
lence, perhaps we are looking for a chimera. 

This has recently been enforced by Hunt [2011]: 
. . . But there are good reasons why the answer to the 
big question that Landau and Batchelor raised about 
whether there is a general theory of turbulence is prob- 
ably 'no'. 

Due to our new theory which draws some profit 
namely from the two-fluid ideas of Landau we are more 
optimistic. But skeptical thoughts remain in view of 
the extraordinary low measurement accuracy of the uni- 
versal constants of turbulence compared with the ex- 
tremely high precisision of the fundamental constants 
of physics like e.g. the vacuum speed of light, or the 
mass of the proton [Fritzsch, 2009]. The challenge is 
possibly the dynamic character of turbulence which is 
in the best case characterized as a stable steady state 
rather than a static property like the electron's electric 
charge. 
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